Dynamics of subduction initiation with different evolutionary pathways

Changes of plate motion may have induced subduction initiation (SI), but the tectonic history of SI is different from one subduction zone to another. Izu-Bonin-Mariana (IBM) SI, accompanied by strong backarc spreading and voluminous eruption of Boninites, contrasts with the Aleutians which shows neither. Using finite element models, we explore visco-elasto-plastic parameters and driving boundary conditions for SI evolution. With an imposed velocity, we find three different evolutionary modes of SI: continuous without backarc spreading, continuous with backarc spreading and a segmented mode. With an increase in the coefficient of friction and a decrease in the rate of plastic weakening, the amount of convergence needed for SI increases from ∼20 to ∼220 km, while the mode changes from segmented to continuous with backarc spreading and eventually to continuous without backarc spreading. If the imposed velocity boundary condition is replaced with an imposed stress, the amount of convergence needed for SI is reduced and backarc spreading does not occur. These geodynamic models provide a basis for understanding the divergent geological pathways of SI. First, IBM evolution is consistent with subduction of an old strong plate with an imposed velocity which founders causing intense backarc spreading and Boninitic volcanism. Second, the New Hebrides SI is in the segmented mode due to its weak plate strength. Third, the Puysegur SI is in the continuous without backarc spreading mode with no associated volcanic activities. Fourth, the Aleutians SI has neither trench rollback nor backarc spreading because the slab is regulated by constant ridge-push forces

Implementation and application of adaptive mesh refinement for thermochemical mantle convection studies

Numerical modeling of mantle convection is challenging. Owing to the multiscale nature of mantle dynamics, high resolution is often required in localized regions, with coarser resolution being sufficient elsewhere. When investigating thermochemical mantle convection, high resolution is required to resolve sharp and often discontinuous boundaries between distinct chemical components. In this paper, we present a 2-D finite element code with adaptive mesh refinement techniques for simulating compressible thermochemical mantle convection. By comparing model predictions with a range of analytical and previously published benchmark solutions, we demonstrate the accuracy of our code. By refining and coarsening the mesh according to certain criteria and dynamically adjusting the number of particles in each element, our code can simulate such problems efficiently, dramatically reducing the computational requirements (in terms of memory and CPU time) when compared to a fixed, uniform mesh simulation. The resolving capabilities of the technique are further highlighted by examining plume‐induced entrainment in a thermochemical mantle convection simulation

On the location of plumes and lateral movement of thermochemical structures with high bulk modulus in the 3-D compressible mantle

The two large low shear velocity provinces (LLSVPs) at the base of the lower mantle are prominent features in all shear wave tomography models. Various lines of evidence suggest that the LLSVPs are thermochemical and are stable on the order of hundreds of million years. Hot spots and large igneous province eruption sites tend to cluster around the edges of LLSVPs. With 3-D global spherical dynamic models, we investigate the location of plumes and lateral movement of chemical structures, which are composed of dense, high bulk modulus material. With reasonable values of bulk modulus and density anomalies, we find that the anomalous material forms dome-like structures with steep edges, which can survive for billions of years before being entrained. We find that more plumes occur near the edges, rather than on top, of the chemical domes. Moreover, plumes near the edges of domes have higher temperatures than those atop the domes. We find that the location of the downwelling region (subduction) controls the direction and speed of the lateral movement of domes. Domes tend to move away from subduction zones. The domes could remain relatively stationary when distant from subduction but would migrate rapidly when a new subduction zone initiates above. Generally, we find that a segment of a dome edge can be stationary for 200 million years, while other segments have rapid lateral movement. In the presence of time-dependent subduction, the computations suggest that maintaining the lateral fixity of the LLSVPs at the core-mantle boundary for longer than hundreds of million years is a challenge

Synthesis And Characterization Of Magnetite And Magnetite-Epoxy Polymers Nanocomposites And Their Thermal And Electrical Behaviors [QC176.8.N35 T161 2007 f rb].

Nanozarah magnetit (Fe3O4) disintesis melalui pemendakan beralkali. Kesan pelbagai aditif dan persekitaran akues-organik terhadap saiz, keseragaman saiz dan morfologi zarah magnetit dikaji. Magnetite (Fe3O4) nanoparticles are synthesized via alkaline precipitation. The effect of different additives and aqueous-organic environment on the size, size uniformity and morphology of magnetite particles is studied

Subduction initiation at relic arcs

Although plate tectonics is well established, how a new subduction zone initiates remains controversial. Based on plate reconstruction and recent ocean drilling within the Izu-Bonin-Mariana, we advance a new geodynamic model of subduction initiation (SI). We argue that the close juxtaposition of the nascent plate boundary with relic oceanic arcs is a key factor localizing initiation of this new subduction zone. The combination of thermal and compositional density contrasts between the overriding relic arc, and the adjacent old Pacific oceanic plate promoted spontaneous SI. We suggest that thermal rejuvenation of the overriding plate just before 50 Ma caused a reduction in overriding plate strength and an increase in the age contrast (hence buoyancy) between the two plates, leading to SI. The computational models map out a framework in which rejuvenated relic arcs are a favorable tectonic environment for promoting subduction initiation, while transform faults and passive margins are not